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We numerically study the wetting (adsorption) transition of a polymer chain on a disordered 
substrate in 1 + 1 dimension. Following the Poland-Scheraga model of DNA denaturation, we 
(N . use a Fixman-Freire scheme for the entropy of loops. This allows us to consider chain lengths 

f*i ' of order N ~ 10 5 to 10 6 , with 10 4 disorder realizations. Our study is based on the statistics of 

loops between two contacts with the substrate, from which we define Binder-like parameters: their 
Fj^ ■ crossings for various sizes N allow a precise determination of the critical temperature, and their 

finite size properties yields a crossover exponent <j> — 1/(2 — a) ~ 0.5. We then analyse at criticality 
the distribution of loop length I in both regimes I ~ O(N) and 1 <C I <C N, as well as the finite- 
size properties of the contact density and energy. Our conclusion is that the critical exponents 
for the thermodynamics are the same as those of the pure case, except for strong logarithmic 
corrections to scaling. The presence of these logarithmic corrections in the thermodynamics is 
related to a disorder-dependent logarithmic singularity that appears in the critical loop distribution 
in the rescaled variable A = l/N as A — > 1. 
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■ defined in equation p|. via the analysis of loop statistics between two contacts with the substrate. 



I. INTRODUCTION 



| The effect of disorder on the wetting transition in dimension 1 + 1 has attracted a lot of interest in the last twenty 
years and remains a rather controversial issue 0,0- The wetting model that we consider here is defined as follows. 
The (impenetrable) substrate is located at z = 0. The polymer chain has N monomers, and the position z a of 
, monomer (a) satisfies z a > 0, with z\ = z/v = 0. The partition function of the model reads 



(RW) 



where H — ^2^ =1 s a ^z a ,o- In equation the sum runs over all random walks (RW) with |z a+1 — z a \ = ±1 and 
C^l ' & ~ TcbT * s ^ ne mverse temperature. The contact energies (e a ) are independent quenched random variables. We 
study here the binary distribution (e = with probability p and e — £q with probability 1 — p). 

On the analytical side, efforts have focused on the small disorder limit [HQ- Since the pure wetting transition has a 
specific heat exponent a pure = 0, the disorder is marginal according to the Harris criterion ||. Based on perturbative 
calculations, Ref finds a marginally irrelevant disorder (i) the quenched and annealed critical temperatures coincide 

S(ii) the quenched critical properties are the same as in the pure (or annealed) case, up to subdominant logarithmic 
corrections. Other studies have concluded that that the disorder is marginally relevant 0,0,01 (i) the quenched and 
• annealed critical temperatures differ by a term which has an essential singularity in the disorder strength [2| (h) the 
critical behavior is governed by some non-trivial disordered fixed point. On the numerical side, the same debate on 
the disorder relevance took place. Numerical studies for flat and exponential disorder distributions or for binary 
' disorder distribution fj have concluded that the critical behavior was indistinguishable from the pure transition. 
^ , On the other hand, the numerical study of for binary disorder pointed towards a negative specific heat exponent 
(a < 0). Finally, the study of Gaussian disorder has been interpreted as an essential singularity in the specific 
heat, that formally corresponds to an exponent a = — oo. This paper aims at clarifying the situation for the problem 



II. POLAND-SCHERAGA MODEL OF THE WETTING TRANSITION 



A. Model and observables 

The wetting model of eq. (|TJ is equivalent to the Poland-Scheraga description of DNA denaturation 0, . The 
relation between the two problems is made apparent if one interprets the coordinate z as the relative coordinate of 
the two DNA strands. We accordingly define a forward partition function Zf(a) for a chain of a monomers, with 
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Zi = Zq = 0. From equation we get 



Q-2 



Z f (a)=e-^^2Z f (a'W(a';a) (2) 



where Af(a': a) is the the partition function of a loop going from a' to a. The asymptotic expansion of A/"(a'; a) is 
given by |10| 

Af(a';a)~a 2 a - a ' f(a-a') (3) 

where Co is a constant and /(x) = is the probability to return to the substrate after x steps. In our model, c = I. 
Other values of c are of interest in the DNA denaturation problem In a similar way, we define a backward 

partition function Zf,(a), defined as the partition function of a chain of N — a monomers, with z a = zn = 0, which 
satisfies 

N 

Z b (a) = e-P e * J2 Z b (a')M(a;a') (4) 

a'=a+2 

In these notations, the partition function Z of equation Q is given by Z — Zf(N) = and the probability for 

monomer a to be adsorbed on the substrate is 

*»> - M f£T < 5 » 

where the factor e^ ea in the numerator avoids double counting of the contact energy at a. The contact density on 
the substrate (a quantity of primary importance in the DNA context) is given by 



1 N 

M T ) = ]v5> (a) (6) 



i=l 



In the pure case, 9n{T) is proportional to the energy. Since this is not true in the disordered case, we also consider 
the contact energy 



1 N 

e N (T) = -J2 £ * P( a ) ( 7 ) 



We will also be interested in Pi 00 p(o, 7), defined as the probability of having a loop between monomers a and 7 on 
the substrate 



p ( n <v\ ^/(q)A^(q; 7)^(7) , s , 

7J 7 /, n ( 8 ) 



B. Numerical implementation 

The above equations, explained in more detail in show that numerical calculations of the partition function 
Z will require a CPU time of order 0(N 2 ). The Fixman-Freire method |l3| reduces this CPU time to 0(N) by 
approximating the probability factor f(x) of equation @ by 

1 I 

f{x) = ^h^ f FF {x) =Y. a * e '" lX ( 9 ) 

i— 1 

In equation @ the number / of couples (a^, 6^) depends on the desired accuracy. The parameters (<Zj, hi) are determined 
by a set of non-linear equations. This procedure has been tested on DNA chains of length up to N = 10 6 base pairs 
[Tj,[i3, and the choice I = 15 gives an accuracy better than 0.3%. We have adopted this value throughout this paper. 

Putting everything together, the model we have numerically studied is defined by recursion equations 121 H for the 
partition functions where the loop partition function J\f(a;a') has been replaced by its asymptotic expression |J3}, 
with the Fixman-Freire approximation (eq. |5J for f(x). 
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III. LOCALIZATION OF THE CRITICAL TEMPERATURE 
A. Loop statistics and Binder-like parameters 

We now define the probability measure M^(l) for the loops existing in a sample of length N as follows: for each 
loop length I, we sum the loop probability Pi 00 p(a, a + I) (eq.JSJ) over all possible origins (a) 

N-l 

M N {l) = Y / Pioo P {a,a + l) (10) 

a=l 

The normalization of this measure over I corresponds to the averaged number of loops in a sample of size N , or 
equivalently to the averaged number N8 N (T) of contacts ®) with the substrate : 

M N = J dlM N (l) =N9 N (11) 

This number is thus extensive Mn oc N in the localized phase (T < T c ), and remains finite as N — * oo in the 
delocalized phase (T > T c ). 

The first moment of the loop measure Mjv(Z) 

_ JdllM N (l) 

< 1 >N = JdlM N (l) (12) 

remains finite as N — > oo in the localized phase (T < T c ), whereas it diverges as < I >n^ N in the delocalized phase 
(T > T c ). We thus introduce the rescaled variable 

and the corresponding probability measure A4n(A) for the loops occupying a finite fraction A = l/N of the whole 
sample. This measure .Mat (A) converges respectively in the thermodynamic limit towards 8(A) in the localized phase 
(T < T c ) and towards 8 (A — 1) in the delocalized phase (T > T c ). At the critical point T = T c , A4n(A) converges for 
large N towards a stable measure Ai c (A) with support A G ]0, 1], in addition to a singular part in <5(A) whose weight 
represents the number of finite loops at criticality. 

To locate the critical temperature, it is thus convenient to introduce the following Binder-like parameters [l6| 

where < I >, < I 2 >, < I 3 > are the first moments of the measure Mjy(Z). In the thermodynamic limit N — > oo, the 
parameter B m (T) jumps from B^T < T c ) = to B^T > T c ) = 1. For finite chain lengths (N 1 ,N 2 , ...), the ratios 
Bnx(T), Bff 2 (T) .. cross at T c at some intermediate value < B(T C ) < 1 given by 

< A > c 

where < A m > c = / dA A m A4 C (A) are the moments of the critical loop measure Ai c (A). The parameter R^(T) has a 
similar behavior with the following crossing value 

R(T C ) - < X ~^0 >C (16) 

Computerwise, the evaluation of the full loop distribution / = 1,..,N via ea . (jlUp requires a time growing as iV 2 . 
To keep a computation time of order O(N), we have chosen to sample the rescaled loop measure Mn(A) on a fixed 
number k max of values A& = k/k max with k = 1,2,.., k max . From this sampling of the loop distribution, we define the 
reduced moments 



fc=i 



(18) 
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from which we build the modified Binder parameters 

< A 2 > fc „„. „fjt_._i <A 3 > fc ,_<A> 



B N V 1 )- ~ I ~ > -"AT U ) - u .2 i 19 ) 

< -A >fc™„„ < A >j 



whose properties are the same as the true Binder parameters (|14|l described above, except for the values of the crossing 
points that now depend on k ma x 

B^»\T C ) = < X * >c ' fc — , fl( fc —)(r c )= <A3 >^ r < X> ^ (20) 
< A >c,fe max < A > c . 



where 



h \ m h 

< A m E ( — ^( A * - — ) ( 21 ) 

k— 1 



We now illustrate these notions on the pure case, before we turn to the analysis of the disordered case. 



B. Loop statistics in the pure case 

In the (±1) random walk model Q with a pure substrate (e a = Eq), criticality corresponds to the condition 
e /3c£o = 2 where the substrate is exactly reflexive Q. The reflexive nature of the substrate at criticality holds more 
generally for pure wetting models from a functional RG analysis [Ll\. This means that the partition function Z/v with 
both ends fixed on the substrate is simply given by the number of random walks returning to the origin after N steps 

Z p N ure (T c ) ~ -J= (22) 



/N 

As a consequence, the critical loop measure ifTUt reads 



i.e. at criticality, there are V7V loops, whose lengths are distributed with the random walk first return probability 

rt0 ~ w (24) 

In terms of the rescaled variable A = l/N, the loop measure (|23|l becomes independent of the size N 

M C (X) * ± (25) 

This means that at criticality, there are a finite number of loops whose length I represents a finite fraction of the size 
N of the sample. This property can be understood as follows : the Levy sum of n independent random variables 
distributed with (|24|1 scales as h + h + ■■■In ~ n 2 : a s a consequence, the number n of loops scales with the size 
N ~ l\ + I2 + ■■■l n of the chain as n~ \/~N. And for Levy sums, it is also well known that the maximal length l max 
among the n terms of the sum N ~ l\ + 1% + ...l n actually represents a finite fraction of the sum 0], i.e. the biggest 
loops indeed occupy a finite fraction of the sample. 

With the measure H25[). the crossing values (|15I16[) are B(T C ) = 1/3 and R{T C ) = 9/5. For our sampling procedure 
with k max terms, the crossing values are given in ea. (|20|l . where the moments (|21|) obtained from (|25|l read 

k m ax / 1 \ ra-3/2 

<A m > cA _=£ -— (26) 

We show in Figure ^ the results of our simulations for the modified Binder parameters i|19[l for k max = 100, 
kmax — 1000, and compare them with the calculations of the full Binder parameters 11141 . which correspond to 
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FIG. 1: 



(a) The Binder parameter 73jy max ' (T) of the pure case for k max = 100 (O)i 1000 (A), (*) and chain lengths up to 



N = 8 • 10 s . The measured crossing values at the common critical temperature T c are in agreement with eqs H20I26H which give 



0.3612.. (fc„ 



100), 0.3413.. (kr, 



1000) 



and | (kr, 



— ). (b) The Binder parameter R 



(T) of the pure case for 



kmax = 100 (O), 1000 (A), f (*) and chain lengths up to TV = 8 • 10 5 

temperature T c are in agreement with eqs I|20I26|I which give 1.6699..(fc ma:E = 100), 1.7584.. {kmax = 1000), and | (k max = ^) 



The measured crossing values at the common critical 

IV' 
2 



kmax — Tf- The crossing temperature T c is indeed independent of k max , and the crossing values are in agreement 
with eqs (|20I26|) . 

Let us now briefly describe the properties of the loop measure off criticality. The finite-size scaling form of the 
partition function Zn with both ends fixed on the substrate is [l^ 

Z^{T)c^=Q({T -T c )Vn) (27) 



where the function Q(x) satisfies 

(i) Q(x = 0) = 1 to recover the critical partition function l|22l) . 

(ii) Q(x — > — oo) = — x e x that corresponds to the localized phase 



Z Pure, T<T) „ rj, _ T)2 N e {T-T c fN (2g) 

In this regime, the loop measure l|10[) reads 

-(t c -t) 2 ; 

Ml <T *{l) ~ N(T C ~T) -^j- (29) 

Af»l/(T-T c ) 2 I 3 / 2 

_ (T c — T) 2 1 

i.e. there exists an extensive number N(T C — T) of finite loops distributed with pi oc (l) = - — 7573 • The Binder 

parameters (|14|l thus converge as N — > 00 towards -Boo(T < T c ) = and Roo{T < T c ) = 3. 
(iii) Q(x — > +00) = l/x 2 that corresponds to delocalized phase 



7 pure ( 

' AT>1/(T-T C )3 (T - T c ) 2 iV 3 / 2 



Zr e (T > T c ) _ , * m ^ - (30) 



The factor I/N 3 / 2 means that the substrate becomes repulsive in the delocalized phase 0, 0|- In this regime, 
the loop measure becomes concentrated on S(X — 1) and the associated Binder parameters 1)14(1 converge towards 
B N ^oo(T > T c ) = 1 and R N ^oo(T > T c ) = 1. 

We show on Figure |3| the loop measures A4^(A) for two sizes (AT = 10 5 and N = 4 • 10 5 ), below, at and above T c . 
At T c , the loop measure is independent of N and corresponds to eq (|25(l . Above T c , an A^-dependent minimum shows 
up. 
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FIG. 2: Log-Log plot of the loop measure Mjy(A) of the pure case for N = 10 5 and N = 4 • 10 5 , and T = T c - 2 (loc), T c and 
T c + 2 (deloc). At T c , the two curves are superimposed. Above T c , a minimum appears. 



C. Disordered case 

We have numerically studied the wetting transition for the following binary distribution for the contact energies 

P(s a ) = (l-p)S(e a -e )+p5(e a ) (31) 

with the three dilution fractions p — 0.25, p — 0.5 and p — 0.75. 

We have used the sampling method explained above for the loop statistics with the factor k max — 1000. We have 
computed the modified Binder parameters B^ max \T) and i?^ max ' (T) (|19J) . for each disorder sample, and for sizes 
N = 10 5 , 2 • 10 5 , 4 ■ 10 5 . Both quantities have been then averaged over 10 4 independent samples (from now on, A 

denotes the disorder average of the quantity A). The crossing of B^ m {T) and R^ max \T) then yield reasonable 
error bars in the localization of T c (see Figure |SJ . 



1. Self- averaging properties 
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FIG. 3: Histogram over 10 4 disordered samples of Mn(1 = : y) at criticality in the case p = 0.5, for sizes N = 10 5 (O), 2- 10 5 
(□), 4- 10 5 (A). 



In disordered systems, extensive quantities are expected to be self-averaging, because spatial averages in a given 
sample are equivalent in the thermodynamic limit to disorder averages. For correlation functions, the situation is 
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more subtle, as discussed in details in J2jj for spin-spin correlations in magnetic systems. Since the loop measure 
Mjv(0 HI U|l is a spatially averaged multiple point correlation function, we have studied its self averaging properties. 
We show in Figure l^the histogram over 10 4 disordered samples of M^(l = ^) a t criticality in the case p = 0.5, for 
sizes N = 10 5 , 2 • 10 5 , 4 • 10 5 : the distribution of Mjy(l = ^) over the samples is more and more peaked around its 
average as N grows. 

Concerning the crossing of Binder parameters, we have also checked that averaging separately the moments in the 
numerator and denominator in eqs. (|14|) gives the same values as the averaged Binder parameters, e.g. 



(32) 



This property was also found for usual Binder parameters in magnetic systems |2l| . 



2. Binder parameters crossings 
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FIG. 4: (a) The averaged Binder parameter B( fe ™^)]v(T), for k max = 1000, (N = 10 5 , 2 ■ 10 5 , 4 ■ 10 s ), and p = 
0.25 (CO) 0.5 (□), 0.75 (A). We have rescaled the contact energy eo so that the different Tc(p) are close. The error bars are 
much smaller than the symbols, (b) The averag ed Binder parameter R( k ™^ N (T), for k max = 1000, (N = 10 s , 2- 10 5 , 4- 10 5 ), 
and p = 0.25 (Q), 0.5 (□), 0.75 (A). 
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FIG. 5: (a) The averaged Binder parameter B$ mal °'(T) , for k max = 1000 and N = 1, 2, 4, 8 • 10 5 and p = 0.5 (b) The averaged 
Binder parameter 7?^ max) (T), for k max = 1000 and N = 1, 2, 4, 8 ■ 10 5 and p = 0.5. 
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In Figure 01 we show the results for the averaged Binder parameters B^ max \T) and R^ max \T) for sizes N = 
1,2,4- 10 5 . The numerical values at the crossing depend on the dilution fraction p (|31() : 

(i) for p = 0.25, we obtain W^)(T C ) = 0.370 ± 0.005 and (T c ) = 1.685 ± 0.005. The critical temperature 
T c = 45.43 ± 0.03 is slightly below the annealed temperature T ann = 45.50 (here £o = —270). 

(ii) for p = 0.5, where we have also studied N = 8 • 10 5 , we obtain £( fc — )(t c ) = 0.430 ± 0.005 and R( k ^(T c ) = 
1.547 ± 0.005, with T c = 45.13 ± 0.03 as compared to the annealed temperature T ann = 45.42. (here £o = —350). 

(iii) for p = 0.75, we obtain B ( fi masc) (T c ) = 0.550 ± 0.005 and i?^ max) (T c ) = 1.360 ± 0.005, with T c = 45.15 ±0.05 
as compared to the annealed temperature T ann = 46.82 (here £o = —515). 

These results show that the crossings of Binder parameters allows to locate precisely the critical temperature. We 
now turn to the analysis of critical properties of various observables. 



IV. STUDY OF CRITICAL PROPERTIES 



A. Distribution of loops of length I ~ O(N) at criticality 




The crossings values (B c , R c ) of the Binder parameters vary continuously with the dilution fraction p (Figure^ • This 
means that the measure M. c {\) of loops occupying a finite fraction A = l/N of the sample at criticality also depends 
continuously on p. We present on Figure these loop measures .M C (A) for p — 0.25, p — 0.5 and p — 0.75 in log-log 
plot, and compare them with the pure case p = 0, which corresponds to the straight line (ln.M c (A) = —(3/2) InA, eq 
In the limit A -> (or InA -> -oo), the measures J\A C (X) become asymptotically parallel to the pure case for 
all p, i.e. 

TOA)^^ (33) 

(see also the more detailed study of finite loops below), but otherwise, the loop measures in the disordered cases are 
qualitatively different from the pure case: a minimum occurs, followed by a weak divergence as A — > 1 (Figure (b)). 
The form of the divergence suggests that it is logarithmic with a p dependent exponent. The simplest form that can 
represent the critical loop measure Ai c (X) on the full interval < A < 1 and that is compatible with all our data is 

«=i( u Fi|F) < 34 » 

The values S(p = 0.25) ~ 0.1, 5(p = 0.5) ~ 0.25 and S(p = 0.75) ~ 0.4 for the exponent <5, and the common value 
C(p — 0.25) ~ C(p — 0.5) ~ C(p = 0.75) ~ 6 for the amplitude C give good fits of (i) the measures M C (X) on 
the whole range A £ [0, 1], with a correct location of the p dependent minimum, and (ii) to the values of the Binder 
parameters crossings shown in Figures 01 and [S] 
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Beyond this numerical evidence, it would be of course interesting to have a theoretical explanation for the appearance 
of this logarithmic singularity in the loop measure near A — > 1 in the disordered case. The only qualitative argument 
we can think of at this stage is the following : in the pure case, we have seen that a minimum appears in the loop 
measure in the delocalized phase (see Figure EJ). In the disordered case, one may argue that the minimum of the 
disordered averaged loop measure M C (X) at criticality comes from the fact that at T c , among the disordered samples 
of size N, a fraction of these samples tend to be slightly delocalized, with a minimum in their loop measure M. c {\), 
whereas the other samples tend to be slightly localized (see Figure[7J|. In other words, if one imagines to associate to 
each sample i a sample-dependent pseudo critical temperature T^(i), as it was done in other disordered systems |22| . 
the presence of the minimum at A m ;„ < 1 reflects the spreading of the pseudo critical temperatures T^(i) around 
the thermodynamic critical temperature T c . The fact that the exponent 8{p) grows with p could be interpreted as a 
consequence of a growing dispersion of the pseudo-critical temperatures TJ? (i) with the strength of the disorder. 
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FIG. 7: At criticality for the case p = 0.5: Log-Log plot of the averaged loop measure M C (X) over 10 4 samples, as compared to 
the loop measure of two particular samples. The upward (resp. downward) oriented curve points towards a delocalized (resp. 
localized) trend 



B. Measure of finite loops at criticality 

In the section above, we have discussed the statistics of loops of length / ~ O(N). We now discuss the statistical 
properties of finite loops, i.e. of length 1 <C / <C N. We have measured the critical loop measure Mjv(Z) for 
N = 1, 2, 4 • 10 5 with fixed values of I ( I = 10 k ; k = 1, 2, 1000). We obtain that the dependence in I is a pure 
power law 

<* j^tj (35) 

with an effective exponent j(N) which increases towards the pure value 3/2 as N increases. For instance for p = 0.5, 
we get 7(7V = 10 5 ) ~ 1.46, -f(N = 2 • 10 5 ) ~ 1.468, -/(N = 4 • 10 5 ) ~ 1.475. This power-law behavior with exponent 
3/2 is moreover in agreement with the asymptotic behavior 133JI for loops of length 0(N) in the limit A — > 0, as it 
should, since the large I behavior of the finite loop measure Mjj c (l) should match the small A — > behavior of O(N) 
loop measure Ad c (X). This requirement actually determines the iV-normalization of finite loops. Indeed, we have 
obtained that at criticality, the measure _M C (A) of O(N) loops is independent of the size N and is well described 
by the form J53J. Via the change of variable I — \N, this leads to the following normalization for the loop critical 
distribution l(33|) 
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C. Contact density at criticality 



As explained at the beginning, the contact density JHJ) is directly related to the normalization of the loop measure 
via l|ll|) . The result i|36|) for the normalization in N of the loop measure thus yields the following finite-size behavior 
for the contact density at criticality 



IN 



(T c )cx 



C(p) 



(37) 



i.e. the leading scaling behavior is the same as in the pure case, but there are strong logarithmic corrections to scaling. 
Plotting (y/~N 0n(T)) for various N thus yields a very poor determination of T c , in marked contrast with its precise 
location through the crossings of the Binder parameters, where these logarithmic corrections are absent. 

We have directly computed 0^{T C ) for N = 10 5 , 2 • 10 5 and 4 • 10 5 , for p = 0.25, 0.5 and 0.75. Our results are 
in agreement with eq. P7|). with the same values of C(p) and S(p) quoted above, just after eq. Q34|). This further 
supports the form 134f) of the critical loop distribution. 



D. Energy at criticality 



We now consider the contact energy J7J. In the binary case, it is closely related to the contact density (jBJl, since 
n == ~\~ f _/y , where is the contact density of diluted sites. At criticality, we expect that r/v scales at most like 
ejq, which implies that the energy eAr(T c ) has the same finite size properties as 9n(T c ) , eq. 1)3 7|l 



e N (T c )oc-=(l 



(In N) s (p) ) 



(38) 



with a coefficient C e {p) < C(p) .The direct measure of the ratios g^|^ c j increases very slowly with N, and typical 
values for N = 2.10 5 are 0.87 (p = 0.25), 0.74 (p = 0.5) and 0.6 (p = 0.75). 



E. Finite-size scaling in the critical region 
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FIG. 8: (a) The Binder parameter B N (T - T c ) of the p = 0.5 case, for N = 10 5 (O), 2 • 10 5 (□), 4 ■ 10 s (A) (b) Master curve 
B I x = (T — T C )\/~N\ of equation 1391 for the same data and symbols. 



We are now interested into the finite-size scaling in the critical region. In terms of the specific heat exponent a, 
the singularity of the free-energy is f(T c ) - f(T) - (T c - T) 2 ~ a . Via hyperscaling ( /(T) - 1/£(T) in dimension 
d = 1), the correlation length along the interface diverges as £(T) ~ 1/(T C — T) 2_Q . According to finite-size scaling 
theory, the appropriate rescaled variable is the ratio N/£(T) between the size N of the system and this correlation 
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length As a consequence, we expect that the Binder parameter B^(T) obtained for various sizes N actually 

only depend on the ratio N/£(T) ~ N(T C - T) 2 ~ a or equivalently 



B N {T) =B(x= (T-T C )JV*) with 



1 



2 - a 



(39) 



We show on Figure[S]the master curve obtained with crossover exponent <fi = 1/2, corresponding to a = 0. 
Considering now the finite size scaling for the energy, equations (|38|) for the energy at T c and the crossover exponent 
found in (|39[) suggest the following form 



eN[T) = 



Go((T c -T)VN) + 



(lniV) 5 (p) 



Gi((T D -T)VN) + 



(40) 



where the ... represent higher order terms. Our conclusion is that the critical exponents of the binary disordered case 
are the same as those of the pure case, except for strong corrections to scaling that are p dependent. 



V. CONCLUSIONS AND PERSPECTIVES 



We have studied the two dimensional wetting transition, in the presence of binary disorder for various dilution 
fraction p. Our analysis is based on the probability measure for the loops of length I existing in a sample of size 
N, with both ends of the chain fixed on the substrate. We have first shown how the introduction of Binder-like 
parameters, built out of the first moments of the loop measure, allows to locate precisely the critical temperature. We 
have then found numerical evidence that the critical loop distribution A4 C (X) in the rescaled variable A = l/N £ [0, 1] 
is not a pure power law (in contrast with the pure case), but contains a logarithmic divergence near A — > 1, with a 
p dependent exponent S(p). Finally, we have explained how this singularity in the loop measure induces very strong 
logarithmic corrections to scaling for the contact density, for the energy, and more generally for thermodynamic 
quantities. 

Our present analysis of the binary disordered case raises the question of the dependence of the critical behavior on 
the disorder distribution. Indeed, we have obtained that the critical loop distribution A4 C (X) varies continuously with 
the dilution fraction p of the binary distribution. More generally, we might expect that M. c (\) could depend upon 
the disorder distribution itself. If this is the case, do the critical exponents differ from those of the binary case? Since 
the example of Gaussian disorder has been interpreted in terms of essential singularities of the Kosterlitz-Thouless 
type, we intend to study other types of disorder distributions to clarify this issue. 

Another interesting direction concerns the effect of disorder when the a priori loop entropy J2J in the Poland- 
Scheraga formulation has an exponent c > 2, in which case the pure transition is first-order [9j. Indeed, in the context 
of the DNA denaturation transition, the binding transition between two pure self-avoiding chains on a cubic lattice 
was found to be first order [2j|. The theoretical explanation that has been proposed is that the self-avoidance 
constraint between denaturated loops and the rest of the chain actually induces an exponent c > 2 for the loop weight 
©. The value c ~ 2.11 has been since measured in Monte-Carlo simulations j24[. In the future, we hope to apply 
our method to the disordered Poland-Scheraga model for the case c > 2, and to compare with the results recently 
obtained by B. Coluzzi 25] via Monte-Carlo simulations of self avoiding walks. 
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